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We study the reliability of extractions of \V US \ based on flavor-breaking hadronic 
t decay sum rules. The "(0,0) spectral weight", proposed previously as a favorable 
candidate for this extraction, is shown to produce results having poor stability with 
respect to sq, the upper limit on the relevant spectral integral, suggesting theoretical 
errors much larger than previously anticipated. We argue that this instability is 
| due to the poor convergence of the integrated D = 2 OPE series. Alternate weight 

choices designed to bring this convergence under better control are shown to produce 
significantly improved stability, and determinations of |V^ S | which are both mutually 
compatible, and consistent, within errors, with values obtained by other methods. 
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, I. BACKGROUND 



Three-family unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix implies 



\V ud \ 2 + \V us \ 2 + \V ub \ 2 = l , (1) 

with the V u }, contribution playing a numerically negligible role 0] . Analyses of Ki e3 incor- 
porating recent updates to the K L lifetime 0, the K + _j3|, K L 0] and K s jH| branching 
fractions, and the K£ 3 form factor slope parameters [6|, together with strong isospin- 
breaking and long distance electromagnetic corrections computed in the framework of 
ChPT @, lead to 

/+(0)|K S | = 0.2173 ±0.0008 , (2) 

which, with the Leutwyler-Roos estimate, /+(0) = 0.961 ± 0.008 [9J] (compatible within 
errors with recent quenched and unquenched lattice results (HI), yields (If 

\V US \ = 0.2261 ±0.0021 . (3) 
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This result is in good agreement with expectations based on unitarit y a nd the most 
recent update of the average of superallowed + — > + nuclear (3 decay [ll| and neutron 
decay 12|| results, \V u d\ = 0.9738 ± 0.0003 (If. The ~ 2a discrepancy observed when 



earlier K decay results were employed thus appears finally to have been resolved. One 
should, however, bear in mind two recent developments relevant to \V u d\'- (i) a new 
measurement of the neutron lifetime, in strong disagreement with the previous world 
averag e |l3l| , and (ii) a Penning trap measurement of the Q value of the superallowed 4e V 



decay |14| in significant disagreement with the average used as input in Ref . , and with 
the potential to raise doubts about current evaluations of structure-dependent isospin- 
breaking corrections [l5|. The potentially unsettled \ V u d\ situation makes alternate (non- 
ius) determinations of \V US \ of interest, both as a means of testing the Standard Model 
(SM) scenario for strangeness-changing interactions, and for reducing errors through 
averaging. Two such alternate methods have been proposed recently. 

In the first, ll^s/Vy is extracted using lattice results for fx/ fn in combination with 
experimental results for r^^/H 71 "/^] |l6| | . With the recently updated MILC rif = 3 



unquenched lattice result, fx/fir = 1-198^;qo6 j 171 ] . the first method yields 

|K S | = 0.22451°;°^ , (4) 

compatible within errors with the determination. 

The second of these proposals involves the analysis of flavor-breaking sum rules em- 
ploying strange and non-strange hadronic r decay data |l8j], and forms the subject of the 
rest of this paper. Existing results, based on the "(0, 0) spectral weight" version of this 
analysis, will be discussed as part of the development below. The discussion to follow 
represents an update and extension of the preliminary results presented in Ref. [Til ]. 



II. V us FROM HADRONIC r DECAY DATA 

With I1^ A .^. the spin J parts of the flavor ij = ud,us vector/axial vector 

correlators, Py/^ the corresponding spectral functions, and Ry/A-ij = r[ r 
v T hadronsvy^ij (7)]/r[r~ — > u T e~u e {^)}, the kinematics of hadronic r decay imply |2 

/■m? nay -i 

Rv/A-,ij = U^V^Sew / — (1 - y T f (1 + 2y T ) p^{s) - 2y T ^ (a) (5) 

Jth m r 1 J 



where y T = s/m 2 T , is the flavor ij CKM matrix element, S E w = 1.0201 ± 0.0003 |2l( 
is a short- distance electroweak correction, and the superscript (0 + 1) denotes the sum 
of J = and J = 1 contributions. Eq. (j3J) is written in such a way that both terms on 
the RHS can be rewritten using the general finite energy sum rule (FESR) relation, 

ds w(s)p(s) = - — ; (f dsw(s)U(s) , (6) 



th 27T2 J\ s \ =so 



valid for any analytic weight w(s) and any correlator II without kinematic singularities. 
Quantities Ryf^y, analogous to Ry/A;ij, are obtained by rescaling the experimental decay 
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distribution with the factor (1 — y T ) k u™ before integrating. The corresponding FESR's 
are referred to as the u (k,m) spectral weight sum rules". Similar FESR's can be written 
down for general weights w(s), for sq < m 2 ., and for the separate correlator combinations 
n vM*ij( s ) and sn yM;ij( s )- The corresponding spectral integrals, f*° ds w(s)p^j A ,. j (s), 
will be denoted Rfj(so) in what follows. In FESR's involving both the J = + 1 and 
J = combinations, the purely J = contribution will be referred to as "longitudinal" . 

With this background, the r-based extraction of V us works schematically as fol- 
lows [lH- Given experimental values for the spectral integrals R^(sq), ij = ud,us, 
corresponding to the same w(s) and same s , the combination 

( s o) = -prr-ja prr^ ( 7 ) 

| v ud\ | ►■us | 

vanishes in the SU(3) flavor limit and hence has an OPE representation, 5Rq PE (so), 
which begins at dimension D = 2. Solving for \ V US \, one has 



IK, 



iRZ(so)/\V ud \ 2 ] - SR% PE (s ) 



At scales ~ 2 — 3 GeV 2 , and for weights used in the literature, the dominant D = 2 term 
in 5Rq PE (so) is much smaller than the leading D = contribution and, as a consequence, 
similarly smaller than the separate ud, us spectral integrals (for physical m s , typically at 
the few to several percent level). The OPE uncertainty, A (5Rq PE (s )), thus produces 
a fractional \V US \ error ~ A (5Rq PE (so)) /2 R™ d (so), much smaller than the fractional 
uncertainty on 5Rop E (so) itself. High accuracy for \V US \ is thus obtainable with only 
modest accuracy for 5Rq PE (sq) provided experimental spectral integral errors can be 
kept under control. 

At present, the absence of a V/A separation of the us spectral data means one must 
work with sum rules based on the observed V+A combination. This combination also re- 
duces the fractional ud spectral integral errors. With present ud spectral data 22[ HI 24|], 



these errors are at the ~ 0.5% level for weights used previously in the literature. The 
much smaller strange branching fraction leads to limited statistics and coarser binning for 
the us spectral distribution 2^, 2(| 23] - The K pole term is very accurately known, but 



errors are ~ 6 — 8% in the K* region and > 20 — 30% above 1 GeV 2 . For wei ghts used in 
the literature, the result is us spectral integrals with ~ 3 — 4% uncertainties |25l 1271 128||. 
Experimental errors on |V^ S | are thus at the ~ 1.5 — 2% level, and dominated by uncer- 
tainties in the us sector. The situation should improve dramatically with the increase in 
statistics and improved K identification available from the B factory experiments. 

A number of points relevant to reducing OPE errors are outlined below. Note that use 
of the V+A sum rules has the added advantage of strongly suppressing duality violation 
at the scales considered ^2£^. Working with weights satisfying w(s = s ) = further 
suppresses such contributions |3(j , as does working at scales so > 2 GeV 2 |3l[ . 

A major, and irreducible, source of OPE uncertainty for "inclusive" sum rules (those 
involving both J = + 1 and J = contributions) is that produced by the bad behavior 
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of the integrated longitudinal D = 2 OPE series. This representation displays badly non- 
convergent behavior, order by order in a s , even at the maximum scale, Sq = m 2 , allowed 
by kinematics (3^. Moreover, for the (k,0) spectral weights, those truncations of this 
series employed in the literature can be shown to strongly violate constraints associated 
with the positivity of the continuum (non-fT-pole) part of Pv+a-us( s ) 01 ■ 

The impossibility of making sensible use of the longitudinal OPE representation ne- 
cessitates working with sum rules based on the J = + 1 combination. Since no complete 
J = 0/1 spin separation of the spectral data exists, a phenomenological subtraction of 
the longitudinal parts of the experimental decay distribution is necessary. This can be 
done with good accuracy because the (very accurately known) tt and K pole terms dom- 
inate the subtraction, for a combination of chiral and kinematic reasons j29l . 34j . Small 



continuum us longitudinal corrections have been constrained phenomenologically, by a 
sum rule analysis of the flavor us pseudoscalar channel |35] and a coupled- channel dis- 
persive analysis of the scalar channel which employs experimental Kir phases, ChPT, 
and short- distance QCD constraints as input j36[. The contribution of the resulting phe- 
nomenological longitudinal continuum us model to the r strangeness branching fraction 
(corresponding to the (0, 0) spectral weight, and s = m 2 ) is < 1% of the total. The 
uncertainty in the bin-by-bin continuum longitudinal subtraction thus plays a minor role, 
even for the (0, 0) spectral weight. Dominance of the longitudinal us continuum by the 
-K"q(1430) and -fT(1460) resonance contributions also means the impact of longitudinal 
subtraction uncertainties on the resulting J = + 1 spectral integrals decreases rapidly 
with decreasing s and is much reduced for weights (like the (k > 0, 0) spectral weights) 
which fall off to zero faster at s = sq than does the (0, 0) spectral weight. 

For the remainder of this paper, we focus on the flavor-breaking combination 

Anoo = ng+« to - 4°£L( S ) . (9) 



A. OPE Input 

The OPE representation of All is known up to dimension D = 6. 
The leading, D = 2, term is given by j3?J 

[An(Q 2 )] D=2 = — ^p- [l + 2.333a + 19.933a 2 + 208.746a 3 + (2378 ± 200)a 4 + • • •] 

(10) 

where a = a s (Q 2 )/ir and m s = m s (Q 2 ), with m s (Q 2 ) and az s (QP) the running strange 
quark mass and coupling in the MS scheme. The 0(a 4 ) term [37( was estimated using 
methods which provided an extremely accurate prediction for the 0(d 3 ) coefficient, and 
similarly reliable predictions of the n^-dependent 0(d 3 m 2 ) coefficients of the electromag- 
netic current correlator all in advance of their explicit calculation. For a and rfu , 
we employ exact solutions corresponding to the 4-loop-truncated (3 and 7 functions j39(, 
with initial conditions for m s .1] and ot s |23| 



m s (2 GeV) = 105 ± 25 MeV, a s (m 2 T ) = 0.334 ± 0.022. 



(11) 
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The forms of the D = 4 and D = 6 contributions are well known, and may be found 
in Ref. j^j. The dominant D = 4 contribution is that proportional to the RG- invariant 
strange quark condensate. This is evaluated using ChPT quark mass ratios [40] , GMOR 
for the light quark condensate, and the conventional estimate r c = (me££) / (m s ss) = 
0.8 ± 0.2 for the ratio of the two condensates. D = 6 contributions are evaluated using 
the vacuum saturation approximation (VSA) and assigned an error of ±500%. D > 6 
terms are assumed negligible. This assumption can be tested for self-consistency since, 
for polynomial weights, w(y) = ^2 m c rn y m (with y = s/s ), integrated D = 2m + 2 OPE 
contributions scale as 1/s™. Neglected, but non-negligible, D > 6 contributions will thus 
show up as unphysical instabilities in the output of a given sum rule (in this case, |V^ S |) 
with respect to sq. Since, when c m ^ 0, the D = 2m + 2 contribution to 5Rq PE (sq) is 
unsuppressed by any additional factors of a s , an so-stability test is particularly important 
for weights, w(y), having D > 6 OPE contributions potentially enhanced through large 
values for one or more of the c m with m > 2. The (2, 0), (3, 0) and (4, 0) J = + 1 spectral 
weights, w^°\y) = l-2y-2y 2 +8y 3 -7y 4 +2y 5 , w^(y) = l-3y+10y 3 -15y 4 +9y 5 -2y 6 , 
and w^'°\y) = 1 — Ay + 3y 2 + 10y 3 — 25y 4 + 24y 5 — lly 6 + 2y 7 , are examples of weights 
having such large higher order coefficients. 

The D = 2 OPE integrals are evaluated using the CIPT prescription |4l| . in which the 
RG-improved expression for An, or its Adler function AD(Q 2 ) = —Q 2 dU(Q 2 )/dQ 2 , is 
used point-by-point along the integration contour. To all orders, the two versions of the 
D = 2 integral are necessarily equal, being related by a partial integration. With both 
All and AD truncated at the same order, however, they differ by terms of higher order 
in a. Our central values employ the 0(a 4 )-truncated RG-improved correlator version (for 
arguments in favor of this choice, see Ref. |37||). 

For the scales employed in this study, both the magnitude and error of 8Rq PE (s ) 
are dominated by the D = 2 contribution. The D = 2 error has two important sources. 
The first is an overall scale uncertainty, associated with the error on the input strange 
quark mass, m s {2 GeV). This uncertainty is ~ 50% for the PDG04 input, Eq. (jllj) . but 
should be reduced considerably by ongoing progress in unquenched lattice simulations. 
The second source of error is the truncation of the D = 2 series. 

The D = 2 truncation error is potentially significant because the series in Eq. ()10j) 
is slowly converging near the spacelike point on the FESR contour. In fact, with high- 
scale determinations of a s {Mz) Q corresponding to an rif = 3 coupling d{Q 2 = m 2 ) ~ 
0.10—0.11, the last three terms in Eq. (jlOj) are actually slowly increasing with order at the 
spacelike point throughout the whole of the kinematically allowed region. Convergence 
of the integrated series will thus typically be slow for weights which emphasize this part 
of the contour. The (k, 0) spectral weights, which involve 2 + k powers of 1 — y, fall 
more and more into this category as k is increased, since, on the contour (y = e ), 
|1 — y\ 2+k oc sin 2+k {9/2) is more and more peaked in the spacelike direction. 

We will use two different monitors of the convergence of the integrated D = 2 series. 
The first involves the difference of the truncated correlator and Adler function versions 
of the series, the second the stability with respect to s of the sum rule output. 

For a series with good convergence, the correlator and Adler function versions of the 
truncated sum should be in good agreement and, moreover, show improved agreement 
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with increasing truncation order. We define r^(so) to be the fractional change in the 
relevant integrated order-/c-truncated D = 2 sum produced by shifting from the correlator 
to corresponding Adler function version. Large values of |r^(so)| and /or an increase of 
l r fc( s o)| with k then signal slow convergence of the D = 2 series [42]. We will take 
twice the sum, in quadrature, of the last included term and the difference between the 
truncated correlator and Adler function versions of the sum as our estimate for the 
D = 2 truncation error. The resulting estimate is considerably more conservative than 
those used previously in the literature. 

Regarding s stability, a truncated series well-converged at s = should remain so 
for some range of s < m^.. Since the exact 6R w (s ) would produce a \V US \ independent 
of so, a well-converged truncation of the integrated D = 2 OPE series should produce 
\V US \ values stable over some interval of sq. If, however, neglected higher order terms 
are actually important at so = rn%, they will be even more so at lower scales, making 
the accuracy of the truncated expression even worse at those scales, and producing an 
unphysical sq dependence to the extracted |V^ S | 43[. The absence of a stability window 



in So thus implies the unreliability of the truncated integrated D = 2 series and/or the 
importance of neglected, but non-negligible, higher dimension terms. 



B. Data Input 

For the ud data we use the ALEPH spectral distribution and covariance matrix |2C 
with overall normalization corrected for the small changes in the e, /i, and strangeness 
branching fractions subsequent to the original ALEPH publication. 

For the us data, we have, unfortunately, been unable to obtain the covariance matrix 
from the OPAL collaboration. The OPAL publication [2?} quotes correlated us spectral 
integral errors only for a range of the (k, m) spectral weights, and only for so = m^- This 
information is insufficient to allow the errors resulting from other choices for either So 
or w(y) to be inferred, precluding implementation of the crucial so-stability test, even 
for the (k, 0) spectral weights. We have thus chosen to work with the somewhat older 
ALEPH data j25|, whose covariance matrix is publicly available. The two data sets differ 
mainly in the values of a small number of the strange branching fractions, a particularly 
important difference being that for t~ — > K~7t + tt~v t . To take into account the changes 



in the branching fraction values, we follow the strategy adopted in Ref. |44j. In this 
approach, the ALEPH distribution for each mode is rescaled by the ratio of the current 
world average to the ALEPH 1999 branching fraction value. The resulting rescaled mode- 
by-mode distributions are then recombined to form the modified total us V+A spectral 
distribution [45]. This scheme should be reliable for modes whose rescalings are close to 
1, but is less clearly so for those, like r~ — > K~7r + 7T~v T , where this is not the case. 

The us K pole spectral integral contribution is fixed by the r[i^ M2 ]-based SM predic- 
tion. This is done because (i) the SM prediction is compatible with the observed r — > Kv r 
branching fraction, but ~ 6 times more precise, and (ii) Eq. (jSJ) already pre-supposes the 
validity of the SM mechanism for hadronic r decays. 

There exist three determinations of the branching fraction B [t~ — > K~7t + tt~v t ], the 
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original ALEPH result, 0.214 ± 0.037 ± 0.029%, and the more recent CLEO and OPAL 
results, 0.384±0.014±0.038% and 0.415±0.059±0.031%, respectively. While the CLEO 
and OPAL results are in good agreement, the agreement with ALEPH is less compelling. 
The OPAL us spectral integral results employed the three-fold average, 0.330 ± 0.028%, 
in setting the overall normalization of the K~7t + 7t~ spectral distribution. Results based 
on this normalization are denoted 'ACO' below. To stress the sensitivity to apparently 
minor changes in the branching fraction values, as well as the importance of improved 
precision, we also present results (denoted 'CO' in what follows) corresponding to the 
alternate rescaling, produced by the average, 0.40%, of the OPAL and CLEO central 
values. The 'ACO'/'CO' branching fraction difference, though only 0.07%, represents 
more than 2% of the ~ 3% total strangeness branching fraction and hence has the 
potential to shift \V US \ by as much as ~ 1% (~ 0.0020). 

III. ANALYSIS AND RESULTS 

In this section we first discuss the existing (0, 0) spectral weight analysis, then present 
alternate determinations, based on weights with improved so-stability for \V US \. 

A. The (0, 0) Spectral Weight Determination of V us 

The (0, 0) spectral weight has been proposed in the literature as a particularly fa- 
vorable case for the \V US \ analysis |l8j]. A key potential advantage is the very close 
cancellation between the weighted ud and us spectral integrals. This manifests itself in 
the OPE representation in suppressed values for the integrated D = 2 OPE series and 
hence similarly suppressed values for the m s -induced D = 2 scale uncertainty. This scale 
uncertainty is the dominant component of the estimated theoretical error in Ref. jlH ]. 
being a factor of ~ 2 larger than the estimated D = 2 truncation uncertainty, and much 
larger than any of the other contribution. The combined theory error, ±0.0009 [3], i s 
swamped by the current ±0.0033 experimental error, but, if reliable, would make the 
(0, 0) analysis a very favorable one for use in determining \V US \ once the much improved 
us spectral data from BABAR and BELLE becomes available. Unfortunately, as we 
will see below, the theoretical uncertainty on 5wQp E {m^.) almost certainly significantly 
exceeds the estimate of Refs. (l8j |. Indeed, we will argue that the convergence of the 
integrated D = 2 (0, 0) spectral weight series is sufficiently bad that the (0, 0) analysis 
is, in fact, an unfavorable one for the extraction of |V^ S |. Fortunately, alternatives exist 
with significantly improved convergence behavior, which allow one to take advantage of 
the general approach proposed in Refs. We return to these in the next subsection, 
after elaborating on the problematic features of the (0,0) analysis. 

The (0, 0)-spectral-weight-based results of Refs. |l8( were obtained using the 0(a 3 )- 
truncated Adler function version of the integrated D = 2 OPE series. The truncation 
error was estimated by combining the magnitude of the last (0(a 3 )) included term in 
quadrature with a measure of the residual scale dependence, the latter obtained by chang- 
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TABLE I: Values of (so), as defined in the text, for various weight choices w(y) and sq = m 2 



Weight 


k = 1 


k = 2 


k = 3 


k = 4 


W J=Q+1 


-0.01 


0.06 


0.20 


0.67 


W\Q 


-0.11 


-0.07 


-0.05 


-0.03 




-0.11 


-0.08 


-0.05 


-0.03 




-0.10 


-0.06 


-0.03 


-0.01 



TABLE II: \V US \ as a function of sq for various FESR weight choices and the ACO and CO 
treatments of the us data, so is given in units of GeV 2 



so 


w ACO 


n^.ACO /...AGO n..AGO 
w 10 w 20 W W 


(0,0) 

w co 


wg° wg° w?° 


2.35 


0.2149 


0.2220 0.2243 0.2201 


0.2172 


0.2236 0.2255 0.2218 


2.55 


0.2167 


0.2218 0.2235 0.2203 


0.2192 


0.2236 0.2250 0.2223 


2.75 


0.2181 


0.2218 0.2230 0.2207 


0.2207 


0.2239 0.2246 0.2229 


2.95 


0.2193 


0.2220 0.2227 0.2211 


0.2219 


0.2243 0.2245 0.2235 


3.15 


0.2202 


0.2223 0.2226 0.2216 


0.2228 


0.2246 0.2246 0.2241 



ing the CIPT scale choice fi 2 = Q 2 to fi 2 = £Q 2 , with 0.75 < £ < 1.5. The resulting 
truncation uncertainty is ^o% °f the 0(a 3 )-truncated Adler function sum. The results of 
Ref . |37] , however, show that the truncated sum is decreased by 47% if evaluated instead 
using the 0(d A ) integrated correlator form. This shift is larger by a factor of nearly 2.5 
than the truncation error estimate of Refs. 0| • 

Further evidence for the slow convergence of the integrated D = 2 (0,0) series is 
provided by the values of r^'°\m 2 ), given in Table |l] The values are not small in 
general, and grow rapidly with increasing k. This increase raises serious doubts about 
any truncation error estimate based on features of the 0(a 4 )-truncated sum, even one, 
like ours, which is more numerically conservative than that of Refs. 0] . 

A final illustration of the unreliability of the convergence of the integrated D = 2, 
(0, 0) spectral weight series is provided by the s -stability results, shown in columns 2 
(ACO case) and 6 (CO case) of Table ITTl No stability window for \V US \ is evident in either 
case, as expected given the indications for poor convergence already discussed above. The 
level of instability is much larger than the total theory error estimated in Refs. 0] , even 
if we restrict our attention to the upper half of the interval displayed in the table. 



B. Alternate Weight Choices 



From Refs. |34j, |37J it is clear that the integrated J = + 1, D = 2 OPE series for 



the (k, 0) spectral weights display rather unfavorable convergence behavior, making such 
weights problematic for use in extracting |V^ S |. Although (k,0) spectral weights occur 
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frequently in treatments of hadronic r decay data, one should bear in mind that one of the 
primary reasons for their introduction, namely the possibility of using them in inclusive 
analyses, is entirely negated by the necessity of avoiding inclusive analyses, which follows 
from the extremely bad behavior of the integrated longitudinal OPE representation. 

In Ref. j23|, the possibility of constructing weights more suitable for use in J = + 1 
non-inclusive sum rules was investigated. These weights were chosen to (i) emphasize 
contributions from regions of the contour showing improved convergence for the D = 2 
All series, (ii) suppress contributions from the region of the spectrum above ~ 1 GeV 2 
where us spectral errors are large, and (iii) control the size of higher order coefficients 
which might otherwise enhance D > 6 contributions. Three such weights, Wiq, Wio, and 
W20 [H|, were constructed, all having profiles on the timelike axis intermediate between 
those of the J = + 1, (0, 0) and (2, 0) spectral weights, and hence similarly intermediate 
relative us spectral integral errors. The much improved D = 2 convergence, compared 
to that of the (k,0) spectral weights, is shown explicitly in Ref. j^ijl. Further evidence 
for this improvement is contained in Table HI which shows good agreement, improving 
with increasing truncation order, between the correlator and Adler function versions of 
the truncated D = 2 sums. The contrast to the (0, 0) spectral weight case is striking. 

As So is decreased, terms in the integrated D = 2 series of 0(a k ), with k > 4, grow in 
size relative to the leading th order term. This growth is least rapid for u> 2 o and most 
rapid for wiq. Though the coefficients multiplying these terms in All are not known, 
this nonetheless indicates that stability for the improved D = 2 convergence will be best 
for W20 and worst for w\o- The r™ 10 (so) values in fact display a cross-over in sign and 
increase in magnitude with increasing k < 4 below sq ~ 2.55 GeV 2 , signalling probable 
deteriorating convergence. We thus base our final results on the highest available scale, 
s = m 2 , and favor the w w and w 2 q analyses over that based on w 10 . 



C. Results 

Much improved sq stability is observed for the weights with improved D = 2 con- 
vergence, particularly w 2 o and Wi . Contributions to the errors on \V US \ for the various 
weights are given in Table HTT1 [ 47| . Sources contributing < 0.0003 theoretical uncertainty 
for all weights considered are not listed explicitly but are included in the total theoretical 
error. Since the results for our favored weights, W20 and wiq, are in excellent agreement, 
and the combined errors are minimized for the latter, we take as our final determination 
the w\o results. Displaying the larger of the asymmetric theory errors, the ACO (CO) us 
data treatments yield the following results, both compatible, within errors, with those of 
Eqs. © and ®: 



0.2223 (0.2246) ± 0.0032 e:C p ± 0.0038 t/l 



(12) 
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TABLE III: Contributions to the error on |V^ S |, at sq — ?7i r , for various FESR weight choices. 



Weight 


us data 


ud data 


m s -scale 5a s 5r c D = 2 truncation 


Theory (total) 


U;(0,0) 
W W 

w w 


±0.0040 
±0.0031 
±0.0028 
±0.0033 


±0.0006 
±0.0006 
±0.0007 
±0.0006 


1:S ±£888? ±0.0000 ±0.0020 
±8:881? ±8:888^ ±0.0005 ±0.0009 
±8:8811 ±:8:888i ±0.0008 ±0.0014 
±8:881! :8:8881 ±0.0004 ±0.0004 


±0.0022 

+0.0038 
-0.0029 
+0.0054 
-0.0041 
+0.0028 
-0.0022 



IV. CONCLUSIONS 

We have shown that the values of \V US \ extracted using the (0, 0) spectral weight sum 
rule display a sizeable instability with respect to So- This instability, combined with 
the results of Ref. |37j . strongly suggests that the true D = 2 truncation uncertainty 
is much larger than previously estimated. Given the level of instability, even our much 
more conservative estimate, shown in row 1 of Table 11111 seems far from being overly 
conservative. We see no plausible way of obtaining a reliable, but more restrictive, 
estimate of this uncertainty. The D = 2 truncation error thus represents a sizeable, and 
irreducible, limitation on the accuracy of the (0, 0) spectral weight determination of \V US \. 

In contrast, for the weights Wio, W20 and wio (whose integrated D = 2 series, by de- 
sign, display improved convergence), good consistency, and much improved So-stability, 
is found. At present the theoretical errors on \V US \ for these weights, shown in Table ITTT1 
are dominated by the m s -scale uncertainty. Near-term improvements in unquenched lat- 
tice simulations should significantly reduce this error. A determination of m s (2 GeV) to 
±5 MeV, for example, would reduce the combined wio theory error to ±0.0013, bringing 
a sub-1% determination of \V US \ easily within reach with the improved B factory data. 
Note that w\o would be favored over W20 because of its smaller truncation error. 

It is also possible to construct weights having, simultaneously, improved D = 2 conver- 
gence and reduced m s -scale sensitivity. We have generated a number of such weights, but 
find they typically weight the region of the spectrum above 1 GeV 2 more strongly than 
do those weights discussed above. As a result, with current us data, they produce very 
large experimental errors, and are not presently useful. We will report on these weights 
elsewhere, once the improved us B factory data has become available and meaningful 
stability tests can be performed. 
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